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Lecture Notes on Fluid Dynamics 

(1.63J/2.21J) 
by Chiang C. Mei, MIT 



1.2 Kinematics of Fluid Motion -the Eulerian picture 

Consider two neighboring stations (not two fluid particles) x and x' at the same instant t, 
where Sx = x' — x is small. The fluid velocity at the two stations are related by 

q(x, t) = q(x, t) + (x - x) ■ Vq{x , t) + 0(x - xf (1-2.1) 

Hence 

Sq(x , t) = q(x', t) - q(x, t) = Sx ■ V(f(£, t) + 0(Sxf (1.2.2) 
Let us introduce the index notation: 

qi=u, q 2 = v, q3 = w; x x =x, x 2 =y, x 3 = z (1.2.3) 

and Einstein's convention: Repeated indices are summed over the range from 1 to 3, and 
the summation symbol is omitted but implied. For example, 
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QiQi = lili = Ql + ll + 93 = Q " Q 

i=l 

Thus we may write (1.2.2) as 

8q i = 5x j ^ : , i = 1,2,3. (1.2.4) 



Now 



% 1 ( dq { t dqA 1 f dq { dq 3 



dxj 2 \dxi dxi I 2 \dxi dx 
Define the rate-of -strain tensor by 



and the vorticity tensor by 

Note that 
and (1-2.4) becomes 

5qi = bxjtij + 5xjQij (1.2.9) 
Let us examine the physics of these terms. 
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1.2.1 Rate-of-strain tensor 



In matrix form, the rate-of -strain tensor is : 
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First, the diagonal terms. It is easy to see that en = du/dx is the rate of stretching per 
unit length in the direction of x, e 2 2 = dv/dy is the rate of stretching per unit length in the 
direction of y, and 633 = dw/dz is the rate of stretching per unit length in the direction of 
z. They are the normal components of the rate of strain tensor. 
Note that 

du dv dw _ . . 

en + e 22 + e 33 = e kk = — + — + — = V • q (1.2.11) 

ox ay dz 

is the rate of volume dilatation due to fluid motion. For a proof, let us consider a cube with 
sides (x, x + Ax), (y, y + Ay) and (z, z + Az). After 5t, the side along x will lengthen from 
Ax to Ax + Ax^St = Ax (l + jfj;5t). Similarly, the side along y will lengthen from Ay to 

Ay (l + §|<5t), and the side along z lengthens from Az to Az (l + §j<5t)- Consequently the 
volume V(t) = AxAyAz will change to 



V(t + 8t) = Ax[l 



du _ \ . 
d-x St ) 




V(t) 
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Hence, the rate of volume dilatation is 

1 V(t + 8t) - V(t) 1 dV 
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Next, the off-diagonal terms. Referring to Figure 1.2.1, consider a plane flow in which |^ 
and Ip do not vanish. In the time interval St the side Ax rotates counterclockwise for an angle 
Mi = ^x 

The total rate of angular deformation is 

50i _ 502 
~6t~~St 



^St. The side Ay rotates counterclockwise for an angle 56 2 
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Figure 1.2.1: Rate of strain tensor components 



Thus ei2 = e xy is a rate of angular deformation, called the rate of shear strain. Other 
components ei 3 and e 23 can be interpreted similarly. 

1.2.2 Vorticity tensor 

The matrix form of Qij is 
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Because of the anti-symmetry, there are only three independent components, which can 
also be used to define the vorticity vector (: 
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Hence 
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Figure 1.2.2: Circulation along a closed circle 



What is the physical meaning of £ ? Consider a plane circular disc A bounded by the 
circle C of radius a, see Figure 1.2.2. By Stokes' theorem 



Now let a — > 0, then, 



or, 



The quantity 
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is the average tangential velocity along the circle. Hence ( n /2 is the average angular speed 
of the fluid circling along C, i.e., the average rate of rotation. The line integral above is also 
known as the circulation. 



